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Abstract: We show that in a large volume scenario of type JIB string or F-theory 
compactifications, single thermal inflation provides only a partial solution to the 
cosmological problem of the light volume modulus. We then clarify the conditions 
for double thermal inflation, being a simple extension of the usual single thermal 
inflation scenario, to solve the cosmological moduli problem in the case of relatively 
light moduli masses. Using a specific example, we demonstrate that double thermal 
inflation can be realized in large volume scenario in a natural manner, and the 
problem of the light volume modulus can be solved for the whole relevant mass 
range. We also find that right amount of baryon asymmetry and dark matter can be 
obtained via a late-time Affleck-Dine mechanism and the decays of the visible sector 
NLSP to flatino LSP. 
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1 Introduction 

One of the most attractive features of low energy supersymmetry (SUSY) is the 
uniflcation of gauge couplings at the scale Mqut ~ 10^^ GeV. The GUT scale is 
two orders of magnitude smaller than the Planck scale, and this hierarchy may be 
understood within a higher dimensional setup of string theory, in which gravity lives 
in a relatively large bulk spacetime, while GUT degrees of freedom are conflned on 
branes [1]. Such scenario then requires a moduli stabilization scheme which would 
stabilize the overall volume modulus at a large value in string unit. 

An attractive scheme realizing large compactiflcation volume is the large volume 
scenario (LVS) of [2], proposed in the context of type IIB string theory. This scheme 
involves two Kahler moduli, Tf, and Tg, where = Iie{Ti,) corresponds to a 4-cycle 
volume modulus determining the bulk Calabi-Yau volume as V ^ t^''^, and Ts = 
Re(Ts) describes a small 4-cycle volume supporting a nonperturbative superpotential 
of the form W^p ~ e~°''^\ The interplay between VFnp and a perturbative correction 
to the Kahler potential suppressed by determines the VEV of Th as 



'b ^ ) 



(1.1) 



which can be exponentially large for a moderately large value of ar^. In LVS, the 
hierarchy between the Planck scale and the GUT scale is given by [3] 



Mgut/Mp ~ 



1/2 



(1.2) 



so Tfe = (9(10^) results in the correct hierarchy between the GUT scale and the Planck 
scale. 

The large VEV of the volume modulus implies that its scalar potential is rela- 
tively flat (at least near the minimum), so the volume modulus is relatively light. 
Specficially, the mass of the large volume modulus in LVS is given by 



which would be in the range of cosmological care unless the gravitino mass 7723/2 is 
heavy enough to ensure > 100 TeV. Cosmologically, huge amount of moduli can 
be produced in the form of a coherent oscillation after primordial inflation, so in 
general light moduli cause disasters due to their too late decays. Depending on their 
life-times, they could destruct light elements formed at BBN, produce too much X 
(7)-rays, distort CMBR or provide too much relic density [4, 5]. One way out of 
this cosmological moduli problem is raising-up of the moduli masses, so that moduli 
can decay before BBN. For this, the volume modulus should be heavier than about 
100 TeV as usual moduli should. For r^, = O(IO^), such heavy volume modulus would 
require m^i2 ^ 10^ GeV, and therefore a careful sequestering of gravity mediated 
SUSY breaking if one wishes to realize low energy SUSY in the visible sector with 
soft SUSY breaking masses msoft ~ 1 TeV [6]. 

On the other hand, considering various possible higher order corrections, either 
in a' or in the string coupling, it is likely to be difficult to realize a sequestering 
which would allow msoft/''^3/2 ~ 10~^ \J^^]- case that an enough sequestering is 
not achieved, the large volume modulus becomes relatively light, and then the most 
compelling solution to the cosmological moduli problem is a late thermal inflation 
[10]. As primordial inflation does, thermal inflation can wipe out pre-existing parti- 
cles, for example, moduli. However it works in a limited sense, since energy density 
of thermal inflation itself becomes a source of moduli reproduction. Depending on 
the mass of moduli, such a reproduction could be fatal, and the large volume scenario 
could be in cosmological trouble. 

In this paper, we study the cosmological moduli problem of LVS for the gravitino 
mass range 




(1.3) 



10^ GeV < 



< 10' GeV, 



(1.4) 



which corresponds to the volume modulus mass range 



10-^ GeV < rur, < 10^ GeV 



(1.5) 



-2- 



for Tf, = (9(10^). We show that, if inflaton decays before the large volume modulus 
starts to oscillate, single thermal inflation can solve the problem only for a limited 
range of rrir^. Particularly, the volume modulus mass range m^-^ = O{10~^ — 1) GeV, 
which would be the case if gravity mediated SUSY breaking is not sequestered well, 
is difficult to be cosmologically viable with single thermal inflation only. On the 
other hand, if inflaton decays just before thermal inflation begins, the abundance of 
moduli can be minimized, so the whole range of rriT-^ can be viable with a single-step 
thermal inflation. However, such a very late decay of inflaton is not typical, so we 
consider a double-step thermal inflation as an alternative solution, and show that 
the problem can be solved in a natural manner. It turns out that correct amounts of 
baryon asymmetry and dark matter can be obtained also via a late-time Affleck-Dine 
leptogenesis [11-15] and the decay of the visible sector NLSP to flatino LSP. 

This paper is organized as follows. In section 2, we briefly review LVS and the 
properties of moduli. In section 3, cosmological moduli problem in LVS is discussed. 
In section 4, it is argued that single thermal inflation is not enough to solve the 
cosmological moduli problem for certain range of the volume modulus mass. We then 
clarify the conditions for a double-step thermal inflation to be realized, while solving 
the cosmological moduli problem for the entire modulus mass range in consideration. 
In section 5, considering a speciflc model for the mediation of SUSY breaking, which 
involves D-term mediation associated with anomalous U{1) gauge symmetry and 
also the conventional gauge mediation, we demonstrate that a double-step thermal 
inflation can be realized in a natural manner and solve the cosmological problem of 
the light volume modulus. We also show that the model can produce correct amounts 
of baryon asymmetry and dark matter too. Section 6 is the conclusion. 



2 Large volume scenario and properties of moduli 

In this section, we review briefly the large volume compactiflcation and properties of 
moduli involved. 

The basic properties of LVS models can be viewed from a simple model in 4D 
supergravity (SUGRA) framework. The model is characterized by the following 
Kahler potential and superpotential [2] 

K = -3\nn + jj-^ + 

W = Wo + Ase-"'^^ (2.1) 

where tj = Re(T/) (J = b, s) are the Kahler moduli determining the size of the 
Calabi-Yau (CY) volume as V = r^^^ — Ts^'^ , and K is expanded in powers of l/u. 

3 /2 3 

The ^-term is the stringy a' correction and given by ^a' = C(3)x(^)/ 2(^5 (27r) 
with x{^) cind Qs being the Euler number of the CY manifold and string coupling. 
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respectively. Wq is the tree-level flux superpotential from the stabilization of the 
dilaton and complex structure moduli. As, the flux dependent constant, can be 
taken to be real value by shift transformation of T^: Tg — )■ + ia. Then, one finds 
that the SUGRA scalar potential at large volume limit is given by 

^ ~ o ■Ml _3 + + o Q/2 \^-^) 



This potential has a SUSY-breaking AdS minimum at 



Since it is known that possible uplift process to get dS vacuum does not change the 
vacuum properties related to our argument, we do not care about that in this paper. 
The masses of the moduli are given by 

, Mp 

m^,, ~ VTia, ~ m^^ ~ m3/2 In , 

"^■3/2 

'^3/2 / 1 

"^'^^-^^MM^r -.-0, m,,=m3/2, (2.4) 

where Mp = 2.4 x 10^^ GeV is the reduced Planck scale, 

0, ^ Mpr57(r,^/^), a, ^ MpIm(T,)/(ryV,^/^), 

0, ~ Mp luTb, a, ~ MpIm(T;,)/(r6) (2.5) 

are the canonically normalized scalar fields of and T^, respectively, ^/^^ and ■j/'h are 
their fermionic superpartners, ^3/2 = e^^^\W\ ~ iW^ol/^b^^ is the gravitino mass. 
A generic feature of moduli masses in LVS is that the mass of the large volume 
modulus is suppressed by a power of I/t^ compared to 7713/2, reflecting the fact that 
the no-scale structure of the scalar potential is approximately preserved in the large 
volume limit. The enhancement or suppression factor, the power of ln(Mp/m3/2), in 
the mass spectrum (2.4) is a consequence of the non-perturbative corrections which 
have a crucial role in moduli stabilization. 

In order to know the real values of the moduli masses, we have to specify the 
vacuum value of Tb and the gravitino mass. In the context of local GUT models, for 
universal tree-level gauge kinetic functions /„, the running gauge coupling constants 

1 /2 

are unified at scale Mqut ~ ^'p/'^b [3) 17]. Therefore, the hierarchy between the 
GUT and Planck scales can be naturally obtained for ~ 10'^ — 10^ as 

1/2 Mqut -.^2 /r, a\ 

rj r.——r. 10. (2.6) 

The gravitino mass is an order parameter of SUSY breaking, and its value is related 
with soft SUSY breaking masses for the visible sector, mgoft = C(lTeV), by SUSY 
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breaking mediation mechanism. In the mixed D-term-gauge mediation which will 
be discussed in our paper, mgoft is just of the order of 7713/2. In other mediation 
mechanism, such as small volume modulus mediation (msoft ~ ^3/2/871"^) [2] whose 
mass hierarchy is similar to that of [16] as a general result for the SUSY breaking 
moduli stabilized by non-perturbative superpotential, and U{1)a threshold correc- 
tions (msoft ~ ^^3/2/ VStt^) [9] in which the soft mass squared is generated at one- loop 
order due to the SUSY breaking mass spectrum of the U{1)a vector superfield, there 
can be hierarchy between 7713/2 and soft terms. On the other hand, as discussed in [6], 
if the visible sector is far from the small 4-cycle governing by r^, and the possible loop 
corrections are quite suppressed, then tree-level gravity mediation induces the soft 
parameters much suppressed compared to the gravitino mass as 777soft ~ ''^3/2 /t^''^- 
* Thus, here we allow the gravitino mass in the range of 

102GeV<m3/2<106TeV, (2.7) 

that gives 

lTeV<77i^, <10^TeV, (2.8) 

10MeV< 771^, < lO^TeV. (2.9) 

The life-time of a modulus is determined by interactions with its superpartners 
and ordinary matter fields localized in small cycles. In the latter case, most important 
parts are the Kahler metric and the gauge kinetic function for the visible sector 

K,is = + ■ ■ ■ , 

fa = ka% + ---. (2.10) 

where $jS are the matter chiral superfields, T„ is the Kahler modulus determining 
the visible sector 4-cycle volume, and (■ ■ ■ ) denotes higher order terms. can not 
be identified as T{,, if not, it provides too small gauge coupling constants. This is the 
phenomenological reason why the matter fields should be localized in a small cycle. 
In this case, the functional form of the Kaher metric is given by Zi ^ yj{Ts,Ty)/Tb. 
Then, it is straightforward to calculate the coefficients of the following interaction 
terms between the large volume modulus and the visible sector fields: 



Ml 



+ (c,M,A"A" + Cim^^iJiij^ + /7.C.) } . (2.11) 



*In this type of model, the typical sfermion mass m^. is of order of m^^ ~ '713/2/'''^^*: while the 
gaugino mass Ma is of O (7713/2 /t"^''^)- In order for rrii^. to be the same order of Ma, we need further 
sequestering as e~^^^Zi ~ 1 + 0(l/r^), where Zi is the matter Kahler metric, or m^. is given by 
idV/dcPb)/Mp + Oirnli^rl) so that (m^ J = ©(m^/^/T^). 
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The visible-sector gauge field strength is denoted by F"'^'^, A" is the superpartner of 
the gauge field, the gaugino. (pi and ipi are the scalar and fermion component of $j, 
respectively, m^. and b^^j are quadratic mass parameters for sfermions coming from 
SUSY breaking and supersymmetric contributions, respectively. Ma is the gaugino 
mass, m^. is the fermion mass for ipi + The model dependent coefficients are 
generically 

1 \ . . ml,_ ^ ml 



= ' ~ 1 ~ 1 or ~ Q ~ 1. (2.12) 

Let us briefly discuss the origin of each terms. At a tree-level, the interactions 
between 0^ and gauge fields are rather suppressed, because fa in Eq. (2.10) does 
not have Th dependence. However, due to the Konishi and the super- Weyl anomaly, 
three point interactions are generated at the one-loop level [17]. Thus Pa is generically 
0{1/Sn^). Ci, Cij and Ca are mostly originated from the soft SUSY breaking mass 
terms that depend on (pb as mgoft ~ Mp/r^'' ~ Mp exp{—nb(pb/Mp), where Ub is the 
constant determined by SUSY breaking and its mediation mechanism, so generically 
they are of order unity. In some cases, there are terms like {dV/d(j)b)/Mp in m^. 
and b^p.. , which are vanishing by equations of motion so that the contribution to the 
soft masses are zero. However, the interation between the modulus and the scalar 
fields are extracted as {d'^V/d'^(j)b) S4>b\^i\'^/Mp ~ m'^^S(f)b\fi\'^/Mp. As a result, Cj, 
and Cij can be enhanced by a factor of m'^^/m'^^f^ for m^^^ ^ 'risoft-^ Q depends on 
the origin of the fermion mass term. Through the Higgs bilinear //-parameter and 
soft masses, the Higgs can interact with (pb- Due to the non-zero vacuum value of 
the Higgs fields, there is a mixing between the Higgs and (pb in the mass eigenbasis. 
Then, the SM fermions can interact with (pb through the modulus-Higgs mixing. 
Higgsinos can directly interact with (pb via fiHuHd term in the superpotential or 
from the Kahler potential by Giudice-Masiero mechanism [18]. In addition, the large 
volume modulus can decay to axion ab, the complex counter part of the modulus, 
with a sizable branching fraction via kinetic term. 

All those interactions are suppressed by Planck scale with additional suppression 
factors depending on channels. The decay rates of (pb to daughter particles are given 

by 



^"-''''^ - 487r Ml ' r,,^2a, - (87r2)2 ' T,,^,a, - ' 



t It can be shown that we can take a field basis such that derivative interactions between (jif, 
and ((/?i, ipi) are absent. Then, ah interactions are obtained from the potential terms. This basis is 
good in the sense that we do not need to care about cancelations between the interactions obtained 
from the kinetic part and potential part. 
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provided the decay channels are kinematically allowed. Note that decays to matter 
fields is at most comparable to the axion channel though it depends on models. Note 
also that the branching fraction of the large volume modulus to photons is suppressed 
by 0(10"^) relative to that of axion channel. 

On the other hand, the strength of interactions between the small volume mod- 
ulus and the matter fields is much more enhanced compared to that for the large 
volume modulus. Because, as in Eq. (2.5), 5ts ~ 6(j)s/Mst. with Mgt — Mp/r^^^ being 
the string scale in large volume scenario, the suppression scale of the interaction 
between (p^ and matters is not Planck scale, but string scale [19]. As a result, in case 
that T„ and Tg have a sizable mixing at the tree-level, the small volume modulus 
mostly decays to the gauge boson pair as 

647r M2' 

where 'Js^A'^A'^ is a numerical coefficient of 0(10) taking the number of allowed 
channels into account. 



3 Cosmological moduli problem in LVS 

The abundance of a modulus is highly constrained, depending on its life-time. The 
small volume modulus decays dominantly to SM particles well before BBN epoch, 
hence it is harmless as long as right amounts of baryon number asymmetry and dark 
matter can be obtained. In case of the large volume modulus, the life-time can be 
written as 

/100TeV\^ , , 

r^, = 0.6 sec Bi^^^a.a, — (3.1) 

where Br^^^ata,, is the branching fraction to large volume axions and always sizable 
for all mass range of (pf For m^^ < 100 TeV, the large volume modulus decays 
around or after BBN epoch. In this listed below, there are various constraints 

depending on the life-time, r,/,^. 

• '^<f>b ~ 0-1 sec: The relic abundance is restricted by the effective number of the 
extra relativistic degrees of freedom, AN^s, to which the primordial abundance 
of ^He, D is highly sensitive. Contribution from moduli decay is 

AN P'^^ ^ - ^ ^21 

Br(0, -> SM) X (?./^7.(Td,^J ^ ' ^ 

where = 7/4 and ^'^(Td,,/,^^) are the relativistic degrees of freedom of a left- 
handed neutrino and standard model particles (including neutrinos), respec- 
tively. The observed value, AiVeff < 1 (la) from BBN [20] and AN^s = 
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0.851°;^^ (la) from BBN+CMB+LSS [21], poses a constraint, BrgM > 0.9. 
Such a sizable branching fraction can be achieved if Cj (cjj) ~ 'm?^^/m1^^^ [22]. 
Otherwise, for Ci{cij) ~ 1, BrsM "C 1, and moduh should be a sub-dominant 
energy component of the universe when they decay. Actually, the constraint 
from the extra radiation contribution to the SM radiation is applicable for the 
moduli decaying in a present Hubble time though has to be replaced to 
X (4/11)^/^ for T^^ > 1 sec. However, the constraint is weaker than others 
described below. 

• 0.1 sec < r<^j^ < 10^^ sec: The direct interaction between decay products of 06 
and the background light elements causes significant changes in the abundances 
of light elements from the predictions by standard BBN [23]. If r^^^ > 10^ sec, 
it also affects the CMB spectrum as a deviation from the Planck distribution 
[24]. The corresponding constraint is, however, rather milder than that from 
the abundance of ^He/D. 

• 10^^ sec < r0j_ < 10^^ sec: In this case, moduli decay after the recombination 
era. The most stringent constraint is coming from the diffused X(7)-ray [5, 25]. 
For the life-time longer than the age of the universe {r^^ > 10^^ sec), X(7)-ray 
produced in the galactic center should be also considered. 

• '^4>b ~ lO^^ s^c- ■'■^ ^his case, the only constraint is from the present dark matter 
rehc density, ^dm^^ ^0.11 [26]. 

All these constraints are shown in Fig. 1 where we see that, for lOMeV < m^^ < 
100 TeV, the abundance of the light modulus is constrained to be less than 

fi^,/i2<O(10-«-10-2) (3.3) 

that corresponds to 

^ < 0(10-1^ - 10-^°) (3.4) 
s 

where n^^ and s are the number density of the light modulus and entropy density, 
respectively. 

The cosmological moduli problem is that the typically expected moduli abun- 
dance is much larger than the upper bounds shown in Fig. 1. In the very early 
universe, when the expansion rate becomes comparable to the mass of a modulus, 
the modulus is expected to be produced enormously in the form of coherent oscil- 
lation caused by vacuum misalignment. In case of the small volume modulus, the 
initial oscillation amplitude is expected to be string scale, and the modulus never 
dominates the universe until it decays well before BBN. Hence it is harmless. On the 
other hand, the large volume modulus has Planckian initial oscillation amplitude , 
and dominates the universe right after it starts to oscillate or at least when it decays. 

■I- The initial misalignment of might be typically super-Planckian. For such a large field value, 
the potential is of an exponential form. Then, in the universe dominated by radiation it will be 
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Tcp, [sec] 

10" 10^^ lO'' 10*^10^ 1 




10"'' ■ 



J^Q-IU l_l I I I I I I I I I I I 

10"^ 10""^ 10"^ 10"^ 0-1 1 10 10^ 10^ lO'* 10^ 

m^JGeV] 

Figure 1. Constraints on 0,^^h'^ before (ph decays are presented for a wide range of 
the modulus mass. Allowed region is below the lines obtained from the present dark 
matter density (blue), diffuse X(7)-ray (red), CMB distortion (green), and light elements 
abundance (dark brown and brown). From the conservative point of view, we assume 
CHuHa = C'(m^^/6), so that the branching ratio for (pb — >■ SM is sizable in the mass range 
of > 2mHiggs- 



We assume that the expansion rate during the primordial inflation is much larger 
than the mass of the small modulus. The light modulus starts to oscillate when 
H ~ fn^^. Denoted a.s Y = n^^/s, the late-time abundance of the light modulus is 
given by 



Ft 



BB,0 



X 




m 9* (^* ) 



-1/4 



1/2 



for Ti > m^, 
for Ti < m^^ 



(3.5) 



quickly attracted to near its vacuum position [27]. Once A^t, < Mp with respect to the vacuum 
position, the potential can be approximated as a quadratic form, and the modulus oscillates with 
an amplitude of order Planck scale and behaves like a matter. Therefore, we can simply assume 
that the initial misalignment of the light modulus is of order Planck scale when it starts to oscillate 
in a radiation dominated universe. 
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where Tr is the reheating temperature of inflaton and defined as 



7r2 



L90^7.(^r) 



-1 



(F/Mp)^/^ (3.6) 



with Tj being its decay rate, and we assumed that the energy density of the modulus 
is equal to that of infiaton for F/ < m^^. If Tj > m^^, the late-time entropy release 
in the decay of heavy moduli, (ps, can dilute the abundance of the light modulus by 
a factor 



^■H^) [^J [jS) [^J f^-^) 

at most. Hence the late-time abundance of the light modulus after the decay of the 
small volume modulus is 

Ybb = 1bb,o^ ^ 0.l7.-.A,A, (yg^) (3-8) 
where we used ~ 10^ and ln{Mp/ 1713/2) — 4vr^. If F/ < m^^, we find 

since the reheating temperature of infiation should be larger than about 1 MeV for 
successful BBN. As shown in Eqs. (3.8) and (3.9), irrespective of the reheating tem- 
perature of infiation, the abundance of the large volume modulus is too much to 
match the observational constraints, Eq. (3.4). Hence, the light modulus causes a 
disaster in a consistent cosmology unless its abundance is somehow diluted enough. 



4 Thermal inflation 

The most compelling solution to the cosmological moduli problem is thermal infiation 
[10]. In this section, we show that single thermal infiation is only a partial solution 
to the moduli problem in LVS for the mass scale of moduli(Eq. (2.9)), though it 
depends on the reheating temperature of the primordial infiation. We then consider 
a double thermal infiation as a complete solution to the problem, and clarify how it 
works. 

4.1 A single thermal inflation 

The fiaton field, denoted as X, controlling thermal infiation has a zero-temperature 
potential, 

V{X) = Vo-mj^\X\'^ + --- (4.1) 
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where (■ ■ ■ ) represent possible higher order term(s) to stabihze Then, denoting 
the VEV of X as Xq and requiring vanishing cosmological constant at true vacuum, 
one finds 

Vo ~ m^X^ (4.2) 



3 \ -1/4 



1 rriv „ /IT 



X 



7x^, T,^{-g4T,)] {TxMpf (4.3) 



where mx is the physical mass of the flaton. Fx is the decay rate of X, 7x is a factor 
determined by the coupling of X to SM particles and not a function of Xq, and Td 
is the decay temperature of X, which is assumed to be dominated by SM particles 
to re-estabilish radiation background for a successful BBN. Thermal inflation takes 
place when Vq dominates the energy density of the universe while the background 
temperature is still larger than the critical temperature Tc ~ mx at which X is 
destabilized from the origin and ends thermal inflation. 

The epoch of flaton domination follows thermal inflation, and the eventual decay 
of flaton reheats the universe, releasing huge amount of entropy. As the result, moduli 
are diluted by a factor 



A 



9*s{Tc 



30^ 



(4.4, 



Ti^Tr 



However there is some amount of moduli reproduction in the following way. During 
thermal inflation, the moduli potential is of the form 

V{<P,) = lml<P,' + "-^Vo + ... (4.5) 




2 



. + +... (4.6) 

ml Mv 



where Cb is assumed as a constant of order unity (c;, = 1/3 has been used throughout 
this paper). Hence moduli is shifted by the amount of 



and reproduced after thermal inflation with the amount of 



Y,,,^U^,..my'^ (4.8) 



§ In the LVS, it is difBcult to obtain such a flaton field {X) at low energy in the closed string 
moduli sector [28]. Instead, X is considered as a open string mode along with the SM matter fields. 
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Therefore, the late-time abundance of moduh is the sum of Eqs. (3.5) and (4.8), that 
is"^ ^ 

Ytot = [1bb,o + >Ti,o] ^ (4.9) 

In the right-hand side of the above equation, the first and second terms as functions 
of Xq provide decreasing and increasing contributions, respectively. Hence the total 
moduli abundance can have a minimum at Xq = X™™. Depending on the decay rate 
of inflaton, it is minimized as follows. 

Tj > m^^i The abundance of moduli is minimized at X^^^ satisfying 

1ti,o = SFbb.o (4.10) 

with the minimum abundance, 

>^tr = 41bb,o^. (4.11) 

Since Ibb fn^l^"^ and X™"^ oc ^^^^ ^ we see that the minimum abundance of the 
light moduli increases as m^^ decreases. 

Tj < m^^: In this case, thermal inflation can take place ii T > when the energy 
density of inflaton is comparable to Vq. Then, as long as the energy density of inflaton 
is subdominant at its decay (i.e., Tx < T/), 1bb,o is minimized for Tr ~ Vq^^- That 
is, 

T/l/4 

Vbb,o > ^ (4.12) 
Then, taking the minimum value of Eq. (4.12), we find 

5yBB,o = 2FTi,o (4.13) 

at Xo = Xq"^^'^ with 



YtT > T^yBBfi-. (4.14) 



7. 

The abundance of moduli increases as m^^^ decreases in this case too. 

In Eqs. (4.11) and (4.14), one crucial factor affecting the abundance of moduli 
is the decay rate which depends on the coupling of X to SM particles. For an 
intermediate scale Xq, there are two leading possibilities. 

• yU-term coupling: 

W D -^T^X^'HuHd (4.15) 



^Whcn thermal inflation begins, the contribution of the small volume modulus to entropy is 
negligible. 
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The energy density of the flaton X soon after thermal inflation is nearly equally 
distributed to the radial and axial modes. If there is a symmetry under which 
X is charge, and the symmetry is broken only spontaneously, then only the 
radial mode is relevant in our argument (though the axial mode is likely to be 
the axion for strong-CP problem of QCD). The decay rate of |X| in this case 
is 

It gives too high decay temperature to solve the problem of the light moduli. 

If X is stabilized by higher order term(s), for example, a self-coupling, then 
radial and axial modes have masses similar to each other. For mx ~ ma_y < 
msoft with rUa^ being the mass of axial component of X, the axial mode decays 
later than radial mode. The decay rate is given by 

where Nc^i is the color factor for a particle i, B is the soft parameter of the 
Higgs bilinear term, the mass of the CP-odd Higgs and m^^ is the mass of 
a particle coupled to down-type Higgs. 

Hadronic coupling: 

W D XX^^ (4.18) 

where \E' and ^ are assumed SM charged fields. The effective coupling of the 
radial and axial modes of X to SM sector is 

g2 1 

Aeff ~ t|-^^ (4.19) 

lOTT^ Xq 

where g is the gauge coupling of \E' and ^. If there is a global symmetry 
under which X is charged, there will be very light Goldstone boson to which 
\X\ decays dominantly. This is dangerous because BBN requires a universe 
dominated by SM-like radiation. Therefore, global symmetries should be badly 
broken. In this case, radial and axial mode have mass scales similar to each 
other, and the decay rate is given by 

where Nii, is the number of the pairs of \E' + The perturbativity of gauge 
coupling restricts N-^ to be less than about 8 if (^) is a fundamental repre- 
sentation of 5[/(5) [29]. 
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Figure 2. Top: Tj > Bottom: Fj = Vq^"^ /^/SAIp maximizing dilution effect from in- 
flation. Left: a fixed flaton mass, mx = 5 GeV. Right: a fixed flaton VEV, Xq = 10^^ GeV. 
The light blue colored region is excluded by observations, and the moduli relic abundance 
^^f^h? after a single thermal inflation was depicted as various (dotted, dashed, solid) lines 
depending on the choice of parameter set. Dark red and blue lines are respectively the 
case of 7x = 1 and 7x = l/Svr^ in Eq. (4.2). Solid lines correspond to the case of minimal 
decay temperature = 5MeV. Dotted lines of the left panels correspond to the minimum 
abundances of (ph given by Eq. (4.11) and (4.14), respectively. 

Using Eqs. (4.17) and (4.20), we show the effect of single thermal inflation on the 
moduli problem in Fig. 2. In the figure, the red and blue lines corresponds to flaton 
decay via //-term and hadronic interactions, respectively. The solid lines are for 
specific choices of parameters shown in the figure, while dashed lines are for minimum 
abundances of moduli obtained by using X™™ for given sets of mass parameters. It is 
clear that the case of hadronic interaction is more efficient in solving moduli problem. 
Such an efficiency is due to the low decay temperature originated from the smallness 
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of 7x. In the upper panels of the figure, where we assumed F/ > m^^, we notice 
that if mx ^100 GeV only m^^ > 1 GeV can be viable whereas, for mx < 5 GeV, 
m^^ < 10~^ can also be saved. Saving m^^ = O(10~^ — 1) GeV is possible only if Td 
is pushed down to a few MeV. Note that for such a low decay temperature of flaton 
there may be no way to generate right amount of baryon number asymmetry though 
dark matter may be obtained from axions with a large coupling constant. However, 
as shown in the lower panel, if Tj is close to Vq^"^ / y/SMp for which infiaton decay 
has maximal effect on the abundance of moduli, single thermal inflation can solve 
the moduli problem for the whole range of m^^ much easily. In addition, there may 
be a room for baryogenesis, such as late-time Affleck-Dine leptogenesis. 

Consequently, we notice that if T[ ~ V^^'^ / y/3Mp the moduli problem of the 
light moduli in LVS can be solved by a single thermal inflation for the whole relevant 
mass range. However such a possibility is the most optimistic case and might not be 
realized in nature. Hence, generically, a single thermal inflation can solve the moduli 
problem of the light moduli in LVS only partially. Particularly, the mass scale of 
moduli associated with TeV scale gravitino mass is not viable in general. 

4.2 Double thermal inflation 

In this section, as a minimal extension of a single thermal inflation to solve the 
moduli problem completely, we consider a double-stage thermal inflation and clarify 
how it works. Here, we assume Tj > m^^ for simplicity 

Compared to a single-stage thermal inflation, the existence of the second stage 
of thermal inflation is useful for diluting moduli produced before the second thermal 
inflation. Hence, as long as we can manage for the second thermal inflation not to 
reproduce too much moduli, the double thermal inflation would be able to solve the 
cosmological problem of the light moduli. Since the extra dilution we need from 
the second thermal inflation may not be very big, a bit high reheating temperature 
after the second thermal inflaiton might be possible, and this may open a room for 
late-time Affleck-Dine mechanism to work. 

Let us denote the two flat directions which trigger the end of each thermal 
inflation as Xi and X2, and the potential energy of flatons as 



where (...) is the vacuum expectation value (VEV) and Xi^o is the VEV of Xi at 
zero temperature true vacuum. We assume 




(4.21) 



Vi:>V2 




II In the case of Tj < m^^, the minimal abundance of moduh is given in the same way as the case 
of the single thermal inflation. 



Then, Vi and V2 may drive the first and second thermal inflation, respectively. In 
order to have such an order of thermal inflation, we require 



Tc,i > Tc,2 (4.23) 

where is the critical temperature when Xi is destabilized from the origin. If 
this condition is not satisfied, X2 would be destabilized at or before the end of the 
first thermal inflation and the second thermal inflation will not take place. Since 
Tc,i ~ TTLXi for order unity coupling to thermal bath with rriXi being the tachyonic 
curvature of the zero temperature potential along Xi around the origin, Eq. (4.23) 
implies 

mx, > mx2 (4.24) 

In addition to the condition Eq. (4.23), the background temperature when the second 
thermal inflation begins should be larger than Tc^2, otherwise X2 would be destabi- 
lized before V2 starts to dominate and the second thermal inflation would not take 
place. Hence another condition for the second thermal inflation is 

T{t2) > Tc,2 (4.25) 

where T(t2) is the temperature of the radiation which couples to X2 at the epoch 
when the second thermal inflation begins. The radiation density after the first ther- 
mal inflation is given by 

Pr = Pr,bg + (4.26) 

where pr,bg is the background radiation and Apr is the radiation contribution from 
the partial decay of Xi. Soon after the end of the first thermal inflation, one finds 

Apr ^ ^ — — Pi (4.27) 

where Fxi^-sm is the partial decay width of Xi to SM particles which couple efficiently 
to X2 directly or indirectly, and r^i^sM ^ H was assumed. This should dominate 
over pr^bg at least before the time t2- Therefore, we find 

pM ~ ^^^^^2 = Br(Xi ^ SM) (^^9.{T,,^)y\!X/' (4.28) 
where T^^i is defined as 

^^?*(Td,i)Ti = 3r?M| (4.29) 
with Fi being the total decay rate of Xi. Hence Eq. (4.25) is translated to 

r,._ > _L^,.(Tfe))-i|- (4.30) 
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For intermediate scale VEVs of Xi and X2, and ^ ''7^X2 ~ ^c,2; this condition is 
easily satisfied. Therefore, double thermal inflation can be realized easily. 
The number of e-foldings is given by 



N, 



e,tot 



In 



TcA 



4, 
+ - In 

3 



c,2 



(4.31) 



where T{ti) is the temperature when the first thermal inflation begins, and the 
first / second term in the right-hand side is the contribution of the first / second thermal 
inflation. The dilution factors due to entropy release of both stages of thermal 
inflation are given by 



Ai 



A. 



gUT{t2))T{t2f Uih) 
V2 



LBr(Xi ^ XX) 



3/4 



C 



T'3 T 



where 



Ci 
C2 



g*s{T{t2)) ( TT 



fi'*5(Td,2) 



-3/4 



3/8 



The late-time total abundance of moduli is 



tot 



(^BB,0 + ^TIl,o) 



Ai 



where 1bb,o is given by Eq. (3.5), and 



Tli 



45 



TI2,0 



1 

A^ 



c,i 4>-^ p 



with ^ = 1, 2. It is minimized at X™", Xf™ satisfying 



Ftii,o 

with the minimum abundance given by 



Fti2,o — Fbb,o 
^1 



-'tot 



22 Ft 



BB,0 



7 AiAa' 



(4.32) 
(4.33) 



(4.34) 
(4.35) 

(4.36) 



(4.37) 



(4.38) 



(4.39) 



If the decay temperature of the second thermal inflation is high enough, there is a 
chance for the late-time Affleck-Dine leptogenesis to provide a right amount of baryon 
number asymmetry. In order for the mechanism to work, the flaton which triggers 
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the end of the last thermal inflation should trigger the mechanism too, otherwise it 
is difficult to obtain a dynamics triggering the mechanism. This implies that Affleck- 
Dine field should be in a symmetry-breaking phase when the last thermal inflation 
ends, that is 

Tc,AD > T,,2 (4.40) 

and X2 should be able to lift up the AD fleld. We may also have a good dark matter 
candidate, for example the fermion superpartner of X2, although QCD-axion might 
be still a good candidate. 

As a brief summary of this section, a double thermal inflation is possible if 
following inequalities are satisfled. 

VI > 1^2 , T,,i>Te,2 , T(t2)>T,,2 (4.41) 



5 A model 

We now consider a concrete example which can realize a double thermal inflation. 
Our model is characterized by the following superpotential. 

+ X^X"^^ + X^Y^^ + -^X^Y (5.1) 

where VFmssm is the MSSM superpotential without /i-term, ('^) and $ ($) are 
SM charged matter flelds, X, Y, Z and S are gauge- singlets under the SM gauge 
groups. An anomalous global U{1)pq symmetry is introduced to provide a flaton 
field {X) responsible for the first thermal infiation and solve the strong CP problem 
[30] simultaneously. Z is introduced to trigger the second thermal infiation. It 
reproduces the Higgs bilinear /i-term by its intermediate scale VEV. We impose an 
appropriate discrete Z4 symmetry as in Table 1. The coupling of the seesaw operator 
is constrained so that the mass of left-handed neutrino is given by 

"-if 

where f„ is the VEV of up-type neutral Higgs field. All A^s are assumed to be of 
order unity with ~ Mqut- 

Since Z has a non-zero vacuum value, domain walls due to the spontaneously 
broken discrete Z4 symmetry appears and becomes problematic if they are stable. 
However the symmetry might not be exact due to either anomaly or higher order 
symmetry-breaking terms originated from gravitational violation of global symme- 
try. Hence the walls are unstable unless the Z4 is a discrete gauge symmetry without 
anomaly. Collapsing walls produces a gravitational wave background, whose energy 
density is constrained by observations. For an intermediate scale VEV of Z, obser- 
vational constraints requires a symmetry-breaking stronger than anomaly [31]. 
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Table 1. charges under the [/(1)pq and discrete Z4 symmetry 



Field 


HuHd 


LHu 




$$ 


X 




z 


52 


f/(l)pQ 








1 


-3 


-1 


3 








Z4 


2 


2 














1 


3 



5.1 Potential of singlet scalars and gauge mediation 

We assume that S couples somehow to SM-charged non-SM field(s) so that it is never 
destabilized. Ignoring the VEV of Higgs fields, the potential of the singlets is given 

by 

V = Vo + V^c,{X,Y) + Vz{Z) (5.3) 

where Vq is for present vanishing cosmological constant, and Vpq and Vz are the scalar 
potential of the canonically normalized PQ breaking fields and Z field, respectively. 
We find 



PQ 



m'x\X\ 



2 I ry\2 



mp\Y\ 



\xr + 



3Mgut 

x\W 



-X^Y + h.c 



4M, 



Z* + c.c. 



GUT 



M2 ' 



(5.4) 
(5.5) 



where Mqut is the apparent GUT scale for the visible sector in LVS. Soft SUSY 



breaking terms 



A^, and Az) are quiet model dependent. We consider a 



natural set-up realized in string theory [32], where the visible sector volume modulus 
is stabilized by D-term potential from the pseudo anomalous U{1)a gauge symmetry. 
In this case, the U (1)pq can be a part of U {1)a and the corresponding SUSY breaking 
D-term (-D^) is obtained as of the order of gravitino mass [8]. Then, D-term mediated 
soft scalar masses are proportional to the U{1)a charge as m^. = q^Da, where qf 
is the same as the charge under the f/(l)pQ, A-parameters and gaugino masses are 
generically suppressed with respect to Da, so that \Ai^\, \Az\ <^ \fnx\, In this 

set-up, U{1)a charges are imposed as my = —3mx > 0. We then find 



IXI 





Mgut 









1/2 



Y 
X 



1 



6^3 



A. 



mx 



(5.6) 



and also the SUSY breaking auxiliary components of X and Y are 



X 





^ 6 




Y 




Ak 



(5.7) 
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Here a particulary interesting feature is that jY is enhanced by vnxjA^ in the 
hmit \Af\ -C \fnx\- In LVS with ?7(1)^, we have 

\mx\ ^ /dI = O (Stt^A,) , (5.8) 



for which 

1 



~ 6/Dl=0(m3/2). (5.9) 



87r2 Y 

The presence of the superpotential tern 

W D A$y$$ (5.10) 

allows a gauge mediation of the SUSY breaking by F^. If there are A'^^ flavors of 
$ + <l> which form the 5 + 5 of SU{5), the gauge threshold contributions just below 
the messenger scale Mmess = \X<s,Y\ to the soft parameters are given by 



9l 



167r2 Y 



M. = -iV,( jL^). (5.11) 



Then the set-up gives 



/V^TMgotV^' M IX VI ^^^^^ r.i9^ 

UPQ ~ |X| ~ , Mmcss ~ \mY\ ~ — ^===wpQ, (5.12) 

\ K ) OVO-L'a 

and the soft scalar masses are 

m^,(MM)~^, m^^(D)~/D;, m^^(GM) ~ ^^v^, (5.13) 

where MM = moduli mediation mainly due to Tg dependence of Zi in Eq. (2.10), and 
of the order of m^j2/2aTs = 771^/2/87:'' [2], D = D-teim contribution induced by the 
U{1)a [8], GM = gauge mediation (5.11). Since A^ ~ ^DaI'^t^' and all dimensionless 
parameters are of order unity, the above results give the following qualitative pattern 
of mass scales: 



M - ''^'^ 



167r2' 

msoft(D) ~ msoft(GM) ~ 87rVoft(MM) ~ mg/s. (5.14) 

We then have enough parameters to make gauge mediation dominates soft terms and 
all squarks and sleptons have positive masses squared. Since Z and S are neutral 
under the t/(l)A, in this case, the dominant contribution to the soft mass-squared of 
Z and S in LVS is from modulus mediation. Hence at UV input scale, we expect 

~ mcj ~ msoft(MM) ~ ^ (5.15) 



- 20 ~ 



At low energy scale, if Ag = 0{1), RG-running of the soft mass-squared of Z becomes 
strong. As the result, potential around the origin along Z direction can have tachy- 
onic instability, and Z develops intermediate scale VEV which reproduce yU.-term 

Since the gravitino mass 7723/2 is around soft mass scale in the visible sector, the 
large volume modulus mass will be around 

m^^ ~ O.OlGeV- IGeV. (5.16) 

Thus, the one-step thermal inflation is not enough in this region. However, we noticed 
that singlet scalars in the model Eq. (5.1) can have a natural hierarchy in their soft 
mass-squared so that 

mx:^mz, Xo > Zo (5.17) 

and hence 

Tc,x > Te,z, Vx,o > Vz,o (5.18) 

where Vz,o = \Vz{Zo)\ and Vx,o = Vq — Vzfi. Therefore, a double-stage thermal 
inflation can be realized. Additionally, a late-time Affleck-Dine leptogenesis may 
have a chance to work. LH^ flat direction is also neutral under U{1)a, so it does not 
have tree-level D-term mass associated with the U{1)a However, the presence 
of gauge mediation effect from Fy provides LHu a soft scale mass-squared which is 
expected to RG-run to a negative value at low energy scale. Around electroweak 
scale, the absolute value of the mass-squared is typically larger than that of the 
singlet Z, and note that gauge mediation effect on LH^ takes place only after X is 
destabilized. Therefore, we expect the following order of destabilization. 

Tc,x > T.^LH^ > T,^z (5.19) 

LHu is eventually lifted up and stabilized around the origin as Z develops interme- 
diate scale VEV reproducing /i-term. This is a perfect circumstance to realize the 
late-time Affleck-Dine leptogenesis [11-15]. 

In the following subsections we will describe how the cosmological problem of the 
light modulus can be solved or ameliorated by the double thermal inflation appearing 
in the model Eq. (5.1). We also check if right amounts of baryon asymmetry and 
dark matter can be obtained for the parameter set solving moduli problem. 

** S could also have a similar instability while Z is held around the origin. However, the RG- 
running of m| caused by the A5 interaction term is slower than that of m^, so Z can develop large 
VEV, providing large mass to S, while S is still at around the origin. As the result, S can not be 
destabilized. 

^'1^ There is a dominant negative 1-loop contribution [9] which drives the soft mass-squared of 
LHu be positive around the electroweak scale. If this is all, LHu flat direction would be stable 
around the origin even without /x-term contribution. 
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5.2 Dilution of moduli by double thermal inflation 



In our model (Eq. (5.1)), all the directions including X, Y, Z and LH^ flat directions, 
but except moduli, are expected to be held around origin at high temperature after 
primordial inflation. The flrst thermal inflation begins when 
temperature 



Vo ~ Vjfo at a 



1/4 



2/3 



-1/12..1/3 
h,0 ^X,0 



(5.20) 



with Hx = V^JqI (-\/3Mp), and ends when temperature drops to T^^x and X is 
destabilized from the origin. Subsequently, LH^ flat direction is destabilized. Both 
of X and LH^ condensations contribute to the radiation density by their partial 
decays. The second thermal inflation can occur if the background temperature is 
higher than T^^z when H ~ Hz = Vz^q I ("\/3^p)- The radiation energy density of 
standard model particles at the time is 



V3Br{X SM)TxMpVz^Q 



V3Br{X SM)TxMpVz^Q 



1 + 
1 + 



Tx-s^SM Pr,x(^c,l) 
8 7LHu I ^LH: 



Bi{X -> SM) V rnx 
where Vx = ^x^sm + ^x^aa is the total decay rate of X with 

Tx-i-SM 
^ X^aa 



(5.21) 



1 \ 2 m\r 



1 m\ 



64n X^ ' 



(5.22) 
(5.23) 



and 7x = 0{1), and 



Q '~lLHu 



(5.24) 



with "^LHu = ^(1) and ttiih^ ~ mgoft/voTT^ being the effective mass for the LH^ 
directional fleld at |£o|, which is originated from the fact that LH^ is stabilized by 
radiative effect rather than tree-level seesaw operator. In the last line of Eq. (5.21) 
we used 



^x^o 



Pv,x{tc,i) ~ m%Xa 



(5.25) 
(5.26) 



For intermediate scale of Xq and Zq with nix ~ (8^^) "^z and rriz ~ Tc^z, one 
flnds pr ^ T^zj hence the second thermal inflation can take place. Note that the 
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contribution from LH^ flat direction to radiation after the flrst thermal inflation is 
larger than that from X for rather mild hierarchy between miHu ^ind mx, so the 
additional dilution from LHu decay can be obtained with a factor 



Soon after the second thermal inflation, the energy density of the universe is 
nearly equally distributed to the radial and axial components of Z. Hence eventual 
reheating temperature is determined by the one which decays later than the other. 
For rriz = 0(10) GeV, axial component decays later with a rate given by Eq. (4.17) 
with niaj^ replaced to rria^ (the mass of axial component of Z). Combined with the 
effect of the partial decay of X, the signiflcant amount of entropy release in this 
decay provide a huge dilution to moduli. 

In Fig. 3, we show the late-time total abundance of moduli after double thermal 
inflation. We use mx = 400 GeV, mz = rUaz = 5 GeV ~ rrix/ (Stt^) and Tc^i = mj 
reflecting the constraints from baryon number and dark matter density studied in 
next subsection. The flgure shows that the whole range of m^^^ can be cosmologically 
viable. Therefore, a double-stage thermal inflation provides a complete solution to 
the cosmo logical moduli problem of the light volume modulus in LVS. 

For the parameters we have used in Fig. 3, we found A^e,tot ~ 17. Although 
such an amount of additional e-foldings does not touch CMB scale, some of inflation 
model is incompatible with this [33]. 

5.3 Baryogenesis 

In our model, LH^ flat-direction is destabilized after the end of the flrst thermal 
inflation, but before the second thermal inflation. This make the late-time Affleck- 
Dine leptogenesis work as Z is destabilized and reproduces /i-term [12]. 

Including the dilution due to entropy release in the eventual decay of Z, the 
resulting baryon asymmetry at present is estimated as [12] 



where nz, riL and uad are number densities of \Z\, lepton asymmetry and AD fleld, 
respectively, rriLHu is the mass scale of LHu, and \io\ is the fleld- value of LHu when 
it is lifted up as Z reaches its VEV. For a small CP-violating phase, 5 <^ 1, the 
conserved lepton asymmetry can be expressed as 



where a ~ 0.1 is the efficiency factor of conserving the generated asymmetry [13, 14], 
and rrig is the mass of the angular mode of the LHu direction when it is lifted up 
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(5.28) 



ul = a6mg\£Q 



(5.29) 
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Figure 3. The total abundance of moduli at decay under the double-stage of thermal 
inflation realized in the model, Eq. (5.1), with constraints given by Fig. 1. Parameters 
used are tan/3 = 10, rrigoft = ITeV, B = O.Smgoft, /U = mA = w-softi n^x = 400 GeV, 
5 GeV, T^x = mx, Tcz = mz, Xq = 3 x 10^° GeV and Zq = 2 x lO'^ GeV. 



mz 



m, 



They are taken to provide correct relic densites of baryon and dark matter (Fig. 4). 



and starts to roll in. We find 



ml ~ 



\2 1/ 12 



m^H 



SGeV^ 



0.1 eV/ VlTeV 



106 GeV 



(5.30) 



where mj, is the mass of a left-handed neutrino, v = 174 GeV is the VEV of neutral 
Higgs and we took tan /3 = 10. Hence 



riL 



a 6 



and 



2 X 10"^° 



me 



10"4 



--In) (n) (^) (^) 



Therefore, the obtained baryon asymmetry may match the observation within the 
uncertainties of involved parameters. 

5.4 Dark matter 

The dark matter candidates in our scenario are QCD-axion and the LSP which is 
flatino, the fermionic superpartner of Z, whose mass is expected to be (9(10) GeV. 
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For fpQ ~ 10^" — 10^^ GeV as the PQ-scale in our scenario, the rehc density of QCD- 
axion is subdominant. On the other hand, abundant flatinos can be produced in 
the decays of the next-hghtest-supersymmetric particle (NLSP) denoted as x which 
might be neutrahno. For the parameter set giving a right amount of baryon asym- 
metry, if kinematically allowed, the decay of Z produces too much flatinos, hence it 
should be forbidden [14] and it is the case since m^^j ~ mz is expected. However, 
production from decays of NLSP can provide a right amount of flatinos to match 
observation. The decay rate of the (neutralino) NLSP to flatino is given by 

r. - 3^7.^ (5.33) 

with 7^ = 0{1). The present abundance of axino dark matter for Tz T^^ = 
is then [14] 

We find that for ~ 200 GeV, ma ~ mz ~ 5 GeV, Zq = 0(10^) GeV with Tz = 
0{1) GeV can match the observed dark matter relic density, flah"^ ^ 0.11. 

The abundance of moduli, baryon number asymmetry and relic density of dark 
matter depend crucially and commonly on mz and Zq, hence in Fig. 4 we show 
contours corresponding to the observed baryon number asymmetry (red) and dark 
matter (blue) with the bound on moduli abundance (the boundary of the colored re- 
gion) in {mz, Zq) plane for mx = 400 GeV. In the figure, the light blue colored region 
is excluded due to moduli over-production. As shown in the figure, right amounts 
of baryon number asymmetry and dark matter can be obtained simultaneously for 
m^, < 0.1 GeV. 



6 Conclusion 

In this paper, we examined the cosmological moduli problem in large volume scenario 
in which the overall volume modulus is stabilized at a value of (9(10^) to generate 
the GUT to the Planck scale ratio Mqut/Mp ~ l/rl'^ ~ 10^^^ 

We found that if the primordial infiaton decay rate F/ > m^-^, single thermal 
inflation can solve the cosmological moduli problem only for a limited range of mT-^. 
Particularly, the mass range m^^ = (9(10~^ — 1) GeV is viable only if the decay tem- 
perature of flaton for thermal inflation is about few MeV, for which baryogenesis is 
difficult to be implemented. On the other hand, if primordial inflatons decay just 
before thermal inflation begins, the whole relevant modulus mass range in consider- 
ation, i.e. m^-f, = C(10~^ — 10^) GeV, can be viable with single thermal inflation, 
while allowing a successful late-time Affleck-Dine leptogenesis. 
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Figure 4. Contours corresponding to the observed baryon number asymmetry (red) and 
dark matter (blue) as functions of mz and Zq for m^^ = 0.08 GeV, = 150 GeV and 
?^L/?^AD = 10~^ rriLHu = 200GeV, 4 = 10^ GeV, g*{Tz) = g*{T^) = 100 with all the 
other parameters being the same as Fig. 3. The light blue colored region is excluded due to 
moduli over-production. In the region above the red/blue line, baryon number asymmetry 
and dark matter at present are smaller than observed ones. 

Since such a late decay of primordial inflatons is not typical, and may not be 
realized in nature, we examined an alternative possibility, double thermal inflation 
as a complete solution of the cosmological moduli problem in large volume scenario. 
Considering a concrete example, we showed that double thermal inflation can be 
realized in large volume scenario in a natural manner, and the cosmological problem 
of the light volume modulus can be solved for the whole relevant mass range. We also 
showed that late-time Affleck-Dine leptogenesis can work after the second thermal 
inflation, and flatino LSPs with a mass of few GeV, which are produced through the 
decays of the visible sector NLSP (e.g. neutralino), can provide a right amount of 
dark matter at present. 
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